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The nearby active galaxy IC310, located in the outskirts of the Perseus cluster of galaxies is a 
bright and variable multi-wavelength emitter from the radio regime up to very high gamma-ray 
energies above 100 GeV. Originally, the nucleus of IC 310 has been classified as a radio galaxy. 
However, studies of the multi-wavelength emission showed several properties similarly to those 
found from blazars as well as radio galaxies. In late 2012, we have organized the first contempo¬ 
raneous multi-wavelength campaign including radio, optical. X-ray and gamma-ray instruments. 
During this campaign an exceptionally bright flare of IC 310 was detected with the MAGIC tele¬ 
scopes in November 2012 reaching an averaged flux level in the night of up to one Crab above 
1 TeV with a hard spectrum over two decades in energy. The intra-night light curve showed a se¬ 
ries of strong outbursts with flux-doubling time scales as fast as a few minutes. The fast variability 
constrains the size of the gamma-ray emission regime to be smaller than 20% of the gravitational 
radius of its central black hole. This challenges the shock acceleration models, commonly used 
to explain gamma-ray radiation from active galaxies. Here, we will present more details on the 
MAGIC data and discuss several possible alternative emission models. 
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1. Introduction 

Relativistic jets are found in many astronomical systems, e.g., as part of the so-called active 
galactic nucleus (AGN) which hosts a supermassive black hole its center, surrounded by an ac¬ 
cretion disk and the jets extending perpendicular to the disk. It is thought that particles in the 
jets can be accelerated to extreme energies and thus, contribute to the ultra-high energy cosmic 
rays with energies of E > 10 18 eV. Current key topics of the research of relativistic jets concern 
the connection between the black hole and the jet, the jet structure, and the jet formation process. 
Furthermore, the emission mechanisms at high energies are still a matter of debate. Hadronic and 
leptonic radiation processes leading to the observed gamma-ray emission are being discussed as 
well as the location of this emission region. Some of these topics can be studied by imaging the 
objects in different frequency ranges. However, even in the radio regime, where the best angular 
resolutions can be achieved, observations are yet capable of resolving the region close to the central 
engine where the jet formation process takes place. The best image of this region so far has been 
obtained for M87 by [[!]] at 230 GHz with a size of 11.0 ± 0.8 gravitational radii. 

An alternative way to study these objects in this respect is to investigate the variability of the very- 
high-energy flux. The size of the emission region is constrained by the finite speed of light and 
causality. Hence, fast variability restricts the size of the region to be smaller than the light crossing 
time. In the gamma-ray energy range the imaging air Cherenkov telescopes provide a large collec¬ 
tion area and hence, a possibility to observe flux variations on time scales of minutes. Such rapid 
events have been measured by H.E.S.S. and MAGIC from PKS 2155—304 and Mrk501 in the last 
decade [|j. ||]. These two objects are classified as blazars. According to the unified scheme, blazars 
are AGN for which the angle between the jet-axis and the line-of-sight is small (6 < 10°). At these 
small angles, large Lorentz factors in the outflow will imply high Doppler factors, which have been 
related to high-energy emission 0. 

The lenticular (SO, z = 0.0189) galaxy IC 310 hosts an AGN with a supermassive black hole with 
a mass of 3 x 10 8 M Q inferred from the M-a relation (see 0). This object has been detected at 
high energies above 30 GeV with FermifLAT [^] as well as at very high energies above 260 GeV 
by MAGIC [g §. In the past, the object was considered to be a head-tail radio galaxy with a large- 
scale (~ 400kpc projected) radio jet [9p. However, Very Long Baseline Interferometry (VLBI) 
observations yielded that the jet shows the same direction in the sky from parsec to kiloparsec 
scalcs|[Tc|, |]. The inner jet appears to be blazar-like with a missing counter jet, presumably due 
to relativistically boosted emission (see also Fig. [l], inset). From the total (de-projected) length of 
the kpc-scale jet and the VLBI observations, we can constrain the angle between the jet-axis and 
the line-of-sight to be 10 < 9 < 20 making IC 310 a borderline object between blazars and radio 
galaxies. Further indications for transitional behavior between a radio galaxy and a blazar was 
found in IC 310 in various energy regimes as summarized in 0. 

2. MAGIC observations and results 

As part of an extensive multi-wavelength campaign in late 2012 to early 2013, IC310 has 
been observed by the MAGIC telescopes. In the night of November 12-13, 2012 a bright flare was 
detected (see Fig. [l]) with a significance of ~ 32 a above 300 GeV 0. 
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Figure 1: Combined image of the significance map measured by MAGIC above 300 GeV on November 12- 
13, 2012 of the Perseus cluster of galaxies and the radio map obtained from EVN observations on October 
29, 2012 at 5.0GHz. The figure is reprinted from |||]. 
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Figure 2: Light curve of the flare on November 12-13, 2012 measured by MAGIC above 300 GeV. The pre¬ 
flare from 23:19 to 23:51 (MJD 56243.972-56243.994) is fitted with a Gaussian function with a standard 
deviation of (9.5 ± 1.9) min. For comparison, the two dashed gray lines indicated the Crab Nebula flux and 
five times the flux of the Crab Nebula. The figure is reprinted from |^|. 
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Figure 3: Zoom of the MAGIC light curve of the flare calculated above 300 GeV into the time range 00:57 
(MJD 56244.04) to 01:40 (MJD 56244.07). The black lines show exponential fits to the light curve to the 
rising and decay phases of the peak. The blue line shows the fit corresponding to the slowest doubling time 
necessary to explain the rising part of the flare at confidence level of 95%. The figure is reprinted from ||5|. 

During this night a mean flux of (6.1 ±0.3) x 10 -11 cm~ 2 s _1 above the same energy was 
measured which is four times higher than the high state flux measured during previous measure¬ 
ments in 2009/2010 [j|]. Figure |] shows the light curve above 300 GeV during the night. Fitting 
the light curve with a constant function reveals a z 2 /d.o.f. of 199/58 corresponding to a probabil¬ 
ity of 2.6 x 10" 17 . The light curve shows several sub-structures indicating fast variability on time 
scales of minutes. The pre-flare can be fitted with a Gaussian function with a standard deviation 
of (9.5 ± 1.9) min. Furthermore, we fitted the light curve with an exponential function in order 
to estimate the flux doubling time. For the sub-structure observed in the time range 00:57 (MJD 
56244.04) to 01:40 (MJD 56244.07), a conservative value (95 % confidence level) for the flux dou¬ 
bling time of 4.88 min could be found for the rising phase (Fig. [|). This value was obtained by 
fitting the rising part at 01:29-01:33 (MJD 56244.062-56244.0652) with many exponential func¬ 
tions, each time fixing the doubling time to a different value and calculating the probability of the 
fit. In the frame of reference of IC 310, this doubling time corresponds to 4.9/(1 ±z) min= 4.8 min. 

The averaged spectrum of the night of the flare can be described with a simple power-law 
function d N /d E = /o x (E/lTcV) 1 with a flux normalization of /o = (17.7 ± 0.9 stat ± 2. l sys t) x 
10 l2 TeV 1 cm 2 s 1 at 1 TeV and a photon spectral index of F = 1.90 ±0.04 stat ±0.15 syst be¬ 
tween 70 GeV and 8.3 TeV. Due to the low redshift of IC310, the corrected spectrum due to the 
extragalactic background light changes by only 20% in the flux normalization and ~ 0.1 in the 
index. The spectral energy distribution from the flare and from previous measurements are shown 
in Fig. [ 4 ]. Compared to earlier observations, no significant spectral variability was observed. Our 
findings are similar to the ones in [||]. More details on the study of the spectral variability can be 
found in [|5|] . 

3. Origin of the gamma-rays 

The standard model to explain of the production of gamma-rays in AGN is the so-called shock- 
in-jet model. In this scenario, particles are accelerated to high energies at shock fronts in the jet via 
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Figure 4: Averaged spectral energy distribution obtained by MAGIC during the flare shown together with 
previous measurements [0. §. Open markers/dashed lines show the measured spectrum, whereas filled 
markes/solid lines indicate the spectrum corrected for the absorption due to the extragalactic background 
light according to the model by [11]. For comparison, the spectral energy distribution from the Crab Nebula 
measured by MAGIC is shown in gray [|I0]. The figure is reprinted from [0], 


Fermi acceleration and radiate, e.g., via the synchrotron-self Compton (SSC) process. Due to the 
magnetic field, the electron-positron plasma produces the observed synchrotron radiation at lower 
frequencies and in case of the SSC model inverse-Compton scattering on the same synchrotron 
photon field leads to the observed high-energy emission. The smallest scale of the shock front 
is determined by the smallest diameter of the jet, i.e., the jet base which is itself assigned by the 
event horizon of the black hole. In case of variability observed from PKS 2155—304 and Mrk501, 
the radius of a spherical emission region is smaller than the event horizon Rq. However, these 
observations are still consistent with the standard model due to the relativistic aberration and a 
large Doppler factor [|4j]. 

In case of IC 310, the variability of 4.8 min constrains the radius of the emission region to be 
~ 20% of the event horizon for the aforementioned mass of the black hole. Due to a rather large 
viewing angle between 10° and 20° and low Doppler factors observed with radio telescopes & 
Doppler boosting can hardly solve this problem. Therefore, this fast variability cannot be easily 
explained with the standard shock-in-jet model. Furthermore, in such a small emission region as in¬ 
ferred from the fast variability, yy pair production would lead to the absorption of the TeV gamma- 
rays making the emission impossible to be observed. The optical depth for 10 TeV photons for 
this process can be calculated to be T rr (10TeV) 300(5/4) 6 (T var /lmin) ^Lsyn/lO^ergs J ) 
assuming a Doppler factor of 8 = 4 [|KJ, a variability time scale of T var = 1 min, and a synchrotron 
luminosity of L syn = 10 42 ergs' '. 

We investigated several models that could explain the sub-horizon variability and anisotropic 
radiation (to avoid pair absorption), e.g., mini-jet structures within the jets |[T3j|. jet-cloud/star in¬ 
teractions [14, 15], and magnetospheric models Jl6[]. As extensively demonstrated in [|0], only 
the last model was able to explain all observed properties. Magnetospheric models are originally 
used to explain the high-energy emission from pulsars. Instead of a pulsar, a particle-starved mag¬ 
netosphere is anchored to the ergosphere of the black hole which is assumed to spin maximally. 
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Figure 5: Magnetospheric model for the origin of the gamma-rays around a maximally rotation black hole 
with an event horizon of Rq (black sphere). The blue shape indicates the ergosphere of the black hole 
in which Poynting flux is generated by the frame-dragging effect. Due to the rotation, a charge-separated 
magnetosphere (red) is induced and vacuum gaps (yellow) are formed at the polar regions. Insight these gaps, 
particles can be accelerated to ultra-relativistic energies due to an electric field parallel to the magnetic field. 
The observed gamma-rays can be explained by inverse-Compton scattering and copious pair production via 
low-energy thermal photons originating from the accreted plasma. The figure is reprinted from ||]. 


Electrons and positrons can be accelerated in vacuum gaps due to electric fields that are parallel to 
magnetic fields. Here, the density of the charge carriers needs to be lower than the Goldreich-Julian 
charge density, otherwise a short-circuit occurs. Because the charge carrier density is proportional 
to the accretion rate of the system, the accretion rate needs to be low. Such a low rate can be as¬ 
sumed for IC 310 as the jet power is rather low [i|]. The fast gamma-ray variability observed from 
IC 310 can be explained by a small gap height. This height is calculated to be h ~ 0.2J?g, assuming 
an accretion rate of m ~ 10~ 4 in units of the Eddington rate and a maximally rotating black hole. 
The gaps can be located at the pole of the magentosphere and can point at different angles with the 
jet axis, e.g., at angles of 10° — 20°. The accelerated particles produce an electromagnetic cascade 
leading to a high multiplicity of the electron-positron pairs. After some time, the gap short-circuits 
and the particles move away with the jet. Then, the gap may reopen and produce new flares. The 
emission is produced by either inverse-Compton scattering on a background photon field, or by 
curvature radiation. According to JT7|], such cascading lead to a rather stable energy spectrum with 
a power-law index of 1.9 as observed from IC 310. 

4. Conclusions 

For first time, variability shorter than the event horizon light crossing time has been observed 
that is not conform with a simple shock-in-jet model. Instead, it is consistent with a magentospheric 
model for high-energy emission in AGN. To finally answer which kind of mechanism is responsible 
for ultra-rapid flux variability events, more observations are needed, e.g., with higher sensitivity as 
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provided by the Cherenkov Telescope Array as this potentially allows to measure even faster flares 
as well as rapid flux variations from other, e.g., fainter AGN. Furthermore observing such events 
with simultaneous multi-wavelength coverage can help to constrain emission models. 
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